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ABSTRACT 

We report low resolution near infrared spectroscopic observations of the eruptive star FU Orionis 
using the Integral Field Spectrograph Project 1640 installed at the Palomar Hale telescope. This work 
focuses on elucidating the nature of the faint source, located 0.5" south of FU Ori, and identified in 
2003 as FU Ori S. We first use our observations in conjunction with published data to demonstrate 
that the two stars are indeed physically associated and form a true binary pair. We then proceed to 
extract J and H band spectro-photometry using the damped LOCI algorithm, a reduction method 
tailored for high contrast science with IFS. This is the first communication reporting the high accuracy 
of this technique, pioneered by the Project 1640 team, on a faint astronomical source. We use our low 
resolution near infrared spectrum in conjunction with 10.2 micron interferometric data to constrain 
the infrared excess of FU Ori S. We then focus on estimating the bulk physical properties of FU 
Ori S. Our models lead to estimates of an object heavily reddened, Ay = 8 — 12, with an effective 
temperature of ~ 4000-6500 K . Finally we put these results in the context of the FU Ori N-S system 
and argue that our analysis provides evidence that FU Ori S might be the more massive component 
of this binary system. 



1. INTRODUCTION 

FU Orionis is the prototype of a small class of rare, 
eruptive young stars. "FUOr" outbursts are generally 
interpreted as episodes of significantly enhanced accre- 
tion in the early stages of disk evolution in young low- 
mass T Tauri stars. In 1936, the apparent magnitude of 
FU Orionis brightened by 6 magnit udes and the object 
has been slowly fading since then (Herbig 1966, 1977; 
" 1 " Petrov||T985| |Kenyon et al.||1988| |(Jlarke et| 



Ko: 



otilov 



2005). Only a few other sources have been identi- 



fied to undergo simila r eruptive phenomena with V1057 



Cyg and V1515 Cyg (Kenyon & Hartmann 1991 Lan- 



dolt 1975, 1977) serving along with FU Orionis as the 



classical examples . 

In the optical, the absorption spectra of FUOr objects 
resemble broadened F-G type supergiants along with 
strong wind/outflow signatures in Balmer and other lines 
(e.g. Call, Na I, K I, Li I); and in the near-infrared they 
have spectral absorption features similar to broadened K- 
M supergiant atmospheres. They generally lack the rich 
emission-line spectrum that is characteristic of rapidly 
accreting young low mass stars. FUOr objects also ex- 
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hibit strong infrared through millimeter excesses, consis- 
tent with the presence of an accretion disk and/or a cir- 
cumstellar envelope of dust. Based on these diagnostics, 
about a dozen additional sources have been identified as 
"FU Ori-like" , with a growing number of them known as 
wide b inary systems (e.g. LI 551 IRS 5, |Rodriguez et aT 



1998 ), RNO lB/C[Ke ^bnetalJ ([l99§, A R6A/B[As^ 



pm~£ Reipurth] ( |20U3p ), including FUDri flWang efaT 



2U04p . For a more exhaustive description of the observa- 



tional properties of FU Ori stars, we direct the reader to- 
wards th e review papers by Hartmann & Kenyon (1996) 



and Reipurth & Aspin (2010 



The accretion disk model for FUOr systems, and FU 
Orionis itself in particular, has been constrained using 
interferometric meas urements and develop ed to a so- 
phisticated level by |Zhu et al.| fl2007[ |2Q1Q[ ). These au- 
thors consider broad- band photometry, spectrophotome- 
try, and high dispersion spectroscopy to derive a model 
size for the outbursting region of the inner disk extend- 
ing outward to 0.5-1.0 AU, with an accretion rate of 
M = 2 x 10 -4 Mq yr~ x onto a central star of mass 
0.371 ^. Interfe rometric measureme nts by Eisner fc Hil- 



lenbr and (20 11]), building on th ose o flMaibet et al.| ( |2005p 
and |lVlillan-Gabet et al. (2006), confirm the rotating ac- 
cretion disk model tor this source. Chandra characteriza- 
tion of FU Orionis unveiled a hot heavily absorbed vari- 
able emission component; the variability is thought to be 
the signature of coronal emission in the c lose vicinity of 
the s tar, seen through the accreting gas (Skinner et al. 
|2010[ ). This interpretation seems consistent with the re- 
cently reported presence of strong magnet ic field in the 
innerm ost region of the accretion disk by |Donati et al. 
( 2005 ) , and reinforces the explanation of the outburst 
which invokes a major increase in the surface brightness 
of an accretion disk because of sudden increase in the 



accretion rate \ Hartmann fc Kenyon| [l996 ) . 
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An alternative to the accretion disk model for FUOrs, 
is the h ypothesis of a fast ro tating G-supergiant photo- 
sphere. Herbig et al. (2003) also showed that the opti- 
cal and infra-red spectroscopic features of fast rotating 
G-supergiants can reproduce FUOr observations if the 
"boxy" line profile shapes can be attributed to core line 
emission rather than t o a rotating disk. In an effort to 
reconcile both models |Kravtsova et al.| ([2007) recently 

el with a puffed 



proposed a modified accretion disk mode. 



inner edge, that is c onsistent with the HST /STIS spec- 
tru m of FU Orionis. jHartmann et al.| ( |20lT] ) however fit 
the |Zhu et ah] ( |2007[ ) disk mode l with Ay = 1-5 m ag to 
the same U V spectrum though Zhu et al. (2007) com- 
ment on the blue excess shortward of 2500-2600 A. Since 
the star/disk source is also embedded in a dust shell, ex- 
tinction plays a crucial role in the spectral characteriza- 
tionof FU O rionis at short wavelengths, and |Kravtsova| 
et al. (2007) suggested using the known companion to 
FU Ononis as a reference to obtain a model independent 
estimate of the optical extinction. 

FU Ori S, a stellar companion locat ed 0.4" to the south 
west of FU Ori, was first discovered by Wang et aLl(|2 004) 
using adaptive optics imaging at K-band. |Wang et al.| 
(2004) discussed, based on statistical arguments using 
on J and K colors, that the discovered object was most 
likely gravita tionally bound to FU Ori. Follow up AO ob- 
servations by Reipurth & Aspin (2004) did not establish 
common proper motion due to the large uncertainties on 
the proper motion of FU O rionis (a result of the large dis - 



p77| )) 
FTTTTri 



tance to FU Ori ~ 450 pc Murdin and Penston 
but provided further evidence of the young age of' " 
S based on infrared excess inferred from JHK'L' colors. 
Spectral line diagnostics in the K-ba nd suggest a late G 
or K s pectral type if the inference of |Reipurth fc A spin 
( |2004| ) is correct that Na I and Ca I are present in their 
spectrum while CO absorption is lackin g. Speckle inter- 
ferometry observations at 0.8 micron by |Kravtsova et aL] 
(120081) find A v = 2.2 mag towards FU Ori S, provided 
the star is of spectral type G9 or later. The tempera- 
ture of the best fitting (unreddened) blackbody to the 
K-band continuum was ~ 2500 K and, based on this, it 
was suggested that FU Ori S exhibits a considerable in- 
frared excess. L-band photometry is also available from 



these authors. Beck & Aspin (2012) recently reported 



high spectral resolution J H and K band (R ~ 3000) 
spectra of FU Ori S, obtained with the Gemini NIFS in- 
tegral field spectrograph. Based on line diagnostics they 
showed that FU Ori S was a highly accreting K5 type 
star and suggested that this object might be the more 
massive component o f the system. N-band interferomet- 
ric measurements by |Monnier et al.| ([2009 ) detected the 
companion, which confirms a substantial in frared excess , 
but no flux density was reported. | Skinner et al. (2010) 
identified a centroid offset towards FU Ori S in the soft 
X-ray component of Chandra observations, and argued 
that the most likely explanation was that FU Ori S was a 
weak soft X ray emitter. All of the evidence is consistent 
with FU Ori S being a young low-mass, perhaps K-type, 
T-Tauri star. 

Precision astrometry on follow up AO observation of 
FU Ori, wit h a baseline longer tha n the time elapsed 



Ori and FU Ori S are gravitationally bound. Obtaining 
a well-constrained SED of this object in these spectral 
regions would allow characterization of FU Ori S inde- 
pendently of FU Ori N. In particular, SED measurements 
covering the spectral region expected to be dominated by 
stellar (as opposed to circumstellar) flux can provide an 
extinction estimate. This measurement may or may not 
be applicable to both components of the binary (FU Ori 
N as well as to FU Ori S) depending on the geometry of 
the circumstellar and circumbinary material. Indepen- 
dent knowledge of the extinction to FU Ori N is necessary 
to validate the findings based on accre tion disk models 
( |Kravtsova et al.||2007| [Zhu et al.||2007[ ) in which the ex- 
tinction is either a derived parameter of the models or 
a required input to them. Understanding the extinction 
therefore provides helpful insights concerning the nature 
of the primary (FU Ori) star as well as the geometry of 
its surroundings, including the companion. 

In an effort to further elucidate the nature of FU Ori S 
by quantifying the spectral type and the reddening, and 
to assist in interpretation of the dust shell encircling the 
FU Orionis N/S system, we observed this source with the 
Project 1640 (P1640) Integral Field Spectrograph (IFS). 
We first summarize our observations, our data reduction 
and establish the binary nature of the FU Orionis sys- 
tem based on our data. Then we reconstruct the 0.8 to 
10 micron SED using previously published observations 
and discuss our photometric points in the context of the 
literature. Based on this reconstructed SED we seek to 
characterize the near infrared emission of FU Ori S. Our 
models lead to estimates of an object heavily reddened, 
Ay =8- 12, with an effective temperature of ~ 4000-6500 
K. We furthermore quantify the amplitude of the infrared 
excess produced by circumstellar dust around FU Ori S, 
using our J and H SED as a ph otometric base l ine fo r 



betwe en the Wang et al. (2004) and Reipurth & Aspin 
(2004) epochs can firmly establish whether or not FU 



interferometric data published by Monnier et al. ((20091), 
reprocessed for the purpose of this paper. We finally put 
these results in perspective in the context of the FU Ori 
N-S system and argue that our analysis provides evidence 
that FU Ori S might be the hotter and therefore more 
massive component. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Observations 

FU Ori was observed at Palomar on March 17 th 2009 
within the - 4" x 4" field of view of the P1640 Inte- 
gral Field Spectrograph. The IFS design prioritizes, for 
a given Field Of View, fine spatial sampling over spec- 
tral resolution for (when compared to other m icrolens 
based ISF such as OSIRIS [McElwain et al.| ( [2007| )). This 
design provides the chromatic information necessary to 
discriminate faint point sources from optical artifacts 
( |Sparks fc Ford|2002[ |Crepp et al.l2010[ ) in high contrast 
observations. Specifically, the 1F!S features a moderate 
spectral resolution (R ~ 45) and high spatial resolu- 
tion (X/D = 45-72 mas). Detail o f the instrument can 
be found in Hin kley et al.| ( |2011a[ ). The data consist 
of twelve exposures of 127 seconds with the primary star 
occulted (behind the coronagraphic mask) and two 2 sec- 
ond exposures with the primary star unobscured (off the 
coronagraphic mask). Each exposure produces a set of 
250 x 250 spectra on the infrared detector, each spectrum 
corresponding to the dispersed image of a micro-lens in 
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the focal plane of the IFS. This raw data is then con- 
verted to a cube of 23 image slices, each slice correspond- 
ing to the image of the focal plane microlens array at a 
given wavelength. The wavelength solution of the spec- 
trograph is derived using a set of calibration fram es ob 
tained prior to and during the observing sequence i Zim- 



merman et al.|20TT ). Determination of the spectrograph 



response in the particular case of FU Ori is detailed in 



2.2 
The 



FU Ori system is seen through an optical 
train consisting of the Palomar Adaptive Optics sys- 
tem (PALAO), the A podised Pupil Lyot Coronagraph 



Soummer et al. 2009b and the PI 640 spectrograph 



Hinkley et al. 20iTa|). ] l'he calibrated data (data cubes 



with speckles that have not been suppressed by post- 
processing), with slices integrated over the J and the H 
bands are shown on the top panel of Fig. [I] We detect 
FU Ori S at a separation of 0.491" ± 0.007 and a posi- 
tion angle of 161.2° ± 1.1 with respect to FU Ori N. The 
contrast sensitivity of the data are to AM# = 6.5 mag 
and AMj = 6.2 mag for companions at 0.5" separation, 
somewhat w orse contrast than t ypical of PI 640 contrast 
curves (e.g. |Crepp et al. ( |2010[ )) due to the early stage 
of the data. 

While the companion FU Ori S can be identified in H 
band, it is as bright as the speckles in J and thus can 
only be distinguished from them using the chromatic 
diversity provided by the IFS. Recent upgrades of the 
Palomar high contrast near-inf rared system, includ ing a 
new Adaptive Optics systems (iDekany et al.||200 6 [) and 
an interferometric wavefront calib ration system jvasisht 
et al. 2010 Wallace et al. 2009), will further improve 
the speckle noise level. Direct extraction of the FU Ori S 
spectrum using the raw slices that were integrated to pro- 
duce the top panel of Fig. [I] leads to over-estimating the 
companion's spectro-photometry because of the speck- 
les' outward motion at the location of the object as the 
wavelength increases. 

In order to mitigate for this effect we sought to im- 
prove the Signal to Noise Ratio (SNR) on FU Ori S using 
the PI 640 sp eckle calibration pipeline presented in |Crepp| 
|et al.| ( |2010[ ). This metho d takes advantage of the chro- 
matic diversity of the IFS ( Sparks & Ford 2002 ) and com- 
bines it with optimal Point Spre ad Function (PSF) sub - 
traction algorithms, e.g. LOCI; ( Lafreniere et al.|[2007 ). 
The processed images are shown in the bottom panels 
of Fig. [I] and demonstrate that our reduction approach 
considerably increases the SNR of FU Ori S. 

2.2. IFS spectral calibration 

Calibrating simultaneously the wavelength solution, 
the atmospheric dispersion and the dispersion intrinsic to 
the instrument with the low spectral resolution (R ~ 45) 
of P1640 is a delicate exercise. Indeed, since the telluric 
sky lines are averaged over each wavelength channel, they 
can not be used as a referen ce to derive the wavelength 
solution ( Cushing et al.|2004| ). We thus first calibrate the 
wavelength solution off-line, using a l aser tunable source 
and f ollow the procedure detailed in Zim merman et aI7| 
(2011). Moreover, the relatively small field of view of 
P1640 prevents us from simultaneously obtaining obser- 
vations of a calibration star in a science exposure and the 
presence of the coronagraph prevents us from using the 
primary star as the calibrator. 




Fig. 1. — Top: band averaged P1640 images of the FU Orionis 
system, FU Ori N is occulted by the coronagraphic mask. Top Left: 
J band. Top Right: H band. FU OriS can be identified in the raw H 
band images, but is as bright as the speckles in J band. Photomet- 
ric estimates based on the raw IFS cube will overestimate the flux 
of FU Ori S because of speckle contamination. Bottom: reduced 
images, the detectability of FU Ori S is significantly enhanced. 



We call the Spectral Response Function (SRF) the 
wavelength dependent relationship between the spectrum 
of an astronomical source and its counterpart seen by the 
P1640 detector. In order to derive an accurate SRF, we 
use the non-coronagraphic images of FU Ori N that were 
acquired right before the coronagraphic observation of 
FU Or i S, combined with the published spectrum of FU 



Ori in Greene & Lada (1996). 



to establish that the 



To do so we first need 
SED of FU Ori N has 



and H 

not varied s ince the 1994 epoch published by |Greene &; 
Lada ( 1996). We first derive three spectral response func- 



tions using three well characterized stars (HD 104860 of 
spectral type F8V, HD87696 / A7V, and iLD109011 / 
K2V) that were observed the same night as FU Orionis. 
We then compare our observed spectrum of FU Ori N, 
respectively normalized by each one of these three re- 
sponse functions. Fig. [2] illustrates this comparison and 
shows excellent agreement for the calibrator within 0.02 
airmass of FU Orionis (HD109011). From this agree- 
ment we conclude that the R ~ 45 P1640 J and H SED 
slopes of FU Ori N hav e not varied since 1994 , and thus 
combine the data from Greene & Lada ( |1996[ ) with our 
non-coronagraphic images to derive the final spectral re- 
sponse function we use to characterize FU Ori S, as il- 
lustrated in the bottom right panel of Fig. [2] . Note that 
the effect of using a calibrator at an airmass that differs 
from the FU Orionis system is most severe in the blue 
and red ends of the SED and in the water absorption 
band between J and H. Since the airmass change occur- 
ring during our coronagraphic observing sequence is 0.02 
(i.e. identical to the difference between the FU Ori and 
the HD109011 exposures, bottom left on Fig.|2| it is very 
difficult to constrain the spectral calibration uncertain- 
ties in these regions of SED. We thus chose to discard 
these points for our analysis of the SED of FU Ori S. 
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Fig. 2. — Determination of the atmospheric and instrumental 
Spectral Response Function. Circles: P1640 spectrum of FU Ori 
N, obtained from non-coronagraphic images normalized using three 
different spectral r esponse functions . Solid Line: IRTF spectrum 
of FU Ori N, from [Greene & Lada| Jl996) . Top: SRF calculated 
using HD 104860, chi-squared difference with the IRTF spectrum 
12.5 Second from top: SRF calculated using HD87696, chi-squared 
difference with the IRTF spectrum 11.2 Second from bottom Left: 
SRF calculated using HD109011, chi-squared difference with the 
IRTF spectrum 4.6, note that most of the mismatch resides in the 
water absorption bands. The P1640 spectrum exhibits excellent 
agreement with the IRTF spectrum when the air-mass of the cali- 
brator star is within 0.02 from FU Ori. We thus c onclude that the 
SEP o f FU Orionis has not varied since the 1994 Greene &; Lada 
( 1996 ) and use the non coronagraphic images of FU Ori N to derive 
the bRF. Bottom : The final SFR is computed using FU Ori N as 
a reference star. 



2.3. Extraction of the spectrum of FU Ori S 

To mitigate the contamination of the spectro- 
photometry of FUOri S by residual speckles, we first 
seek to calibrate these residual quasi-static optical ar- 
tifacts using an aggressive PSF subtraction algorithm, 
the Locally Optimized Combination of Im ages approach 
(Crepp et al. 2Q1Q| Lafreniere et al. 2007). For a given 
location in the image, LOCI creates a synthetic reference 
PSF, a weighted sum of images within a collection of ref- 
erence frames, based on a least-square fit of the neighbor- 
ing speckles. In order to preserve as much flux from the 
companion as possible, the least squares fit is calculated 
using a large "optimization" region of the image, with 
area expressed as Na PSF cores, while the actual sub- 
traction is carried out in a smaller "subtraction" zone. 



As discussed in Pueyo et al. (2011), direct photometry 
on images processed using classical LOCI implementa- 
tions, such as the reduced images in the bottom panel 
of Fig. [T] for instance, can produce biases in the spectro- 
photometry. We quantified this effect, by extracting a 
"zeroth-order" spectrum of FU Ori S from the reduced 
data cubes of Fig. [T] and used this spectrum to inject 
a set of synthetic companions in our dataset. The top 
panel of Fig. [3] shows, for a typical azimuthal orienta- 
tion, the injected synthetic spectrum and a series of ex- 
tracted spectra after LOCI reduction. Each extracted 
spectrum was estimated using a reduction with a given 
area of least-squares minimization, Na- Using such a re- 
duction strategy leads to underestimating the flux of the 
companion by a factor of two to three. Moreov er this 
bias i s wave length dependent. As presented in |Pueyo] 
et al.| ( |2011[ ) this flux depletion is a combination of a 
least squares bias common to all LOCI implementations 
(that can be identified as a "grey" gain that does not 
alter neither spectral features nor the SED slope), and a 
spectral cross talk term, specific to integral field spectro- 
graphs. 

We circumvented this proble m by using the "dam ped 
LOCI" approach introduced in Pueyo et al. (2011). d- 
LOCI relies on a least-squares approach similar to the 
one described above but adds a supplemental penalty 
term, that scales with the flux of the discovered compan- 
ion, in the underlying quadratic cost function. This pre- 
serves the flux from faint sources and leads to un-biased 
SED estimated even in the case of a companion buried 
under speckles. The results with synthetic companions 
are shown in the bottom panel of Fig. [3] Injected and 
extracted spectra agree very well for a wide range of d- 
LOCI parameters. We further explored potential reduc- 
tion biases by varying other algorithm parameters that 
scale with the cross-spectral channels contamination, d- 
LOCI reductions yield consistent results for a wide range 
of parameters, with a bias smaller than 5 percent. We 
thus conclude that our final spectrum is unbiased at the 
~ 5 percent level. 

In Figure [4] we present our final spectrum for FU Ori 
S and, for comparison, that of FU Ori N. We derived 
the error bars as the root mean square of the sum of 
three separate terms: the error on the spectral response 
function (e.g. photometric scatter in non coronagraphic 
images), the scatter of bias estimated on synthetic com- 
panions over the ensemble of LOCI parameters explored 
(smaller than 5 percent), and the photometric scatter in 
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Fig. 3. — Quantification of LOCI and d-LOCI bias using syn- 
thetic companions. The spectra are shown before calibration by 
the Spectral response Function. In both panels the dashed line 
corresponds to the synthetic companion injected in the corona- 
graphic PSFs. This "zeroth order" spectrum was obtained using 
aperture photometry on preliminary LOCI reduced images, and 
renormalized to match the band averaged photometry of the raw 
data. The top panel shows the synthetic spectrum extracted af- 
ter LOCI: while the images exhibit the high SNR illustrated on 
Fig. [l] the spectro-photometric signal clearly exhibits an algorith- 
mic flux depletion. Moreover this bias depends upon the reduction 
parameters. We alleviate this problem using a d-LOCI approach, 
bottom panel, which only exhibits a small bias over a wide range 
of reduction parameters. 

the reduced images. The spectrum of FU Ori S is rising 
to the red throughout the J and H bands and appears to 
be featureless. 

3. PHYSICAL ASSOCIATION OF THE TWO 
COMPONENTS OF THE FU ORIONIS SYSTEM 

3.1. Proper Motion of FU Orionis 

The proper motion of FU Or i N is reported both in 
the Carlsberg Meridian Catalog |Evans et al~| ( |2Q0l[ ) FJ 
as pmRA = 0.1167 year ±0.104, pmDEC = (XU8647 
year ±0.0115, and Proper Motion Extended Catalog 
(PPMXL), as pmRA = 0.0144"/ year ±0.0108 and 
pmDEC = 0.0721"/ year ±0.0108 / . Although the proper 

9 http://archive.ast.cam.ac.uk/camc/ 



I s 5 



n 



{mi"' i 



Fig. 4. — J and H band Spectral Energy Distribution of the 
FU Orionis system. Blue Circles: P1640 spectrum of the North- 
ern component. Red Squares: P1640 spectro-photometry of the 
Southern component. 



motion estimates in right ascension have large percentage 
error, the proper motion in declination is highly signif- 
icant. If we adopt a distance of 450 pc to FU Orionis, 
these two estimates yield surprisingly large transverse 
velocities, 183 km/s and 139 km/s respectively. Trac- 
ing back the trajectory of FU Orionis to ~ 1 M yrs us- 
ing these values places this object in the T-Tauri star 
forming M78. While there exists a multitude of obser- 
vational evidence for the young age of FU Orionis, its 
location, seemingly isolated from any star forming clus- 
ter, can be seen as somewhat puzzling. This large proper 
motion could explain this paradox, and in this scenario 
FU Ori could be a runaway star from M78. However the 
large ejection velocity required to explain this scenario 
can only be the result of an unlikely very close encounter 
early in the history of FU Orionis. Moreover such trans- 
verse velocities do not seem comp atible with the radia l 
velocity of FU Ori N measured by Herbig et al. (2003), 
which is consistent with the molecular cloud velocity. Al- 
ternatively these proper motion estimates, based on op- 
tical images, could be biased by the extended optical 
nebulae surrounding FU Orionis. We further explored 
this issue and derived proper motion using astrometric 
estimates based solely on infrared observations, since at 
longer wavelength FU Ori is a true point source. To do 
so, we take advantage of the ~ 10 years baseline between 
2MASS and WISE observations. Using the 2MASS point 
source catalog and the WISE preliminary catalog, we ob- 
tain pmRA = 0.033"/ year ±0.14, pmDEC = 0.025"/ 
year ±0.053. While the mean value of this proper mo- 
tion estimate yields a smaller transverse velocity of 56 
km/s, the large declination uncertainty in the prelimi- 
nary WISE point source catalog, due to a systematic bias 
in WISE astrometric estimates that will be corrected in a 
future release, prevents us from making a definitive con- 
clusion regarding the actual true proper motion of FU 
Orionis. For the purpose of demonstrating that FU Ori 
North and South are co-moving, we will use the small- 
est of these three proper motion values, the 2MASS to 
WISE infrared to infrared estimate. 
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As the o ne year time baselin e betw een th e epochs pre- 



sented in Reipurth & Aspin (2004) and Wang et al. 
(2004) was not sufficiently large to establish or rule 
out physical association between the FU Ori North and 
South, we explore this question using the 2009 epoch 
from PI 640. In our data we detect FU Ori S at a separa- 
tion of 0.491" ± 0.007 and a position angle of 161.2° ± 1.1 
with respect to FU Ori N. The astrometric pupil plane 
grid, necessary to carry out astrometry on occulted coro- 
nagraphic images (Sivara makrishnan et al. ]|2006 Hink 



ley et al. 2010[ |2011b|), was not in the optical path 



during the observing sequence of FU Orionis. Astromet- 
ric estimate were thus carried out using the un-occulted 
reddest images, where the companion is visible, albeit 
at very low SNR. Each of the five images chosen in 
the H band was analysed with the program FITSTARS; 
which uses an iterative blind-deconvolution that fits the 
location of delta functions and their relative intensity 
to the data. The FITSTARS program was presented in 



ten Brummelaar et al. (1996 2000). After throwing out 
any data that failed to converage to a physical solution, 
the position angle and separation was computed by a 
weighted average of the individual values. The weights 
were set equal to the inverse of the RMS residual of 
the fit; a standard output of the FITSTARS program. 
The error bars for position angle and separation were set 
equal to the weighted standard deviation of the results. 

To increase confidence in our astrometric point, we in- 
cluded in our common proper motion analysis a second 
P1640 epoch, obtained in 2011 with the astrometric cal- 
ibration grid in the optical train. We also complemented 
our observation using the 2005 epoch ob tained as part 
of the interferometric survey presented in |Monnier et al.| 
(2009), and re-processed with particular care to extract 
the flux from FU Ori S. The result of our analysis is sum- 
marized in Figure [5j which includes points previously 
published, where both the uncertainties in the relative 
position of both components and the proper motion of 
FUOri N have been combined as an uncertainty in the 
position of FU Ori S. The configuration of this system 
is somewhat cumbersome as the largest source of un- 
certainty is the lack of a well constrained proper motion 
estimate for FU Ori N. Even in the presence of such large 
uncertainty Figure [5] shows that the combination the five 
epochs clearly establishes that FU Ori N and FU Ori S 
are co-moving. We thus conclude that these two objects 
are physically associated. The relative position of the 
binary pair across the six epochs in Figure [5] is given in 
Table [I] Under the assumption of a face on orbit and a 
0.4 Mq mass for FU Ori N, the largest orbital motion de- 
tectable over the 9 year baseline provided by these points 
is 26 micro-arcseconds in the NS direction assuming that 
FU Ori S is massless and 45 micro-arcseconds in the NS 
direction if the mass of FU Ori S is twice the mass of FU 
Ori N. In any case, the astrometric uncertainties in Table 
[I] are much larger than these optimistic projected relative 
motion, and we conclude that the scatter in these values 
is not representative of any significant orbital motion. 



4. RECONSTRUCTION OF THE 0.8-10 /xM 
SPECTRAL ENERGY DISTRIBUTION OF FU 
ORI S 



Uncertainty in the proper motion of FU Orionis + relative uncertainty 
Uncertainty in the relative location of FU Ori S wrt FU Ori N 




0.0 0.2 

RA (arcsec) 



Fig. 5. — Proper motion of the FU Orionis system. The blue 
dots represent the location of FU Ori N across the five epochs 
considered. The star motion was calculated using the astrometric 
difference of the WISE and 2MASs epochs. The uncertainties on 
this motion, have been cumulated with the uncertainty in the rel- 
ative position of FU Ori S with respect to FU Ori N to estimate 
the regions of confidence of the on-sky location of FU Ori S (or- 
ange ellipses). For clarity, we arbitrarily used the 2005 epoch as 
the origin of time for the calculation of the uncertainties in proper 
motion and the 2002 epoch as the zero point of our diagram. Since 
the regions of confidence between the 2002 and the 2003 epochs 
overlap, common proper motion could not be established at the 
time. However the data points introduced in this paper provide 
much greater temporal leverage and demonstrate that FU Ori N 
and FU Ori S are physically associated. 

4.1. Consistency with published J and H band 
observations of FU Ori S 

We first focused on comparing the P1640 band av- 
eraged photometry with the estimates obtained by 
Reipurth fc Aspin (2004). By integrating the data 
points in Fig. Hover the J and H bandpasses, we found 
delta magnitudes between FU Ori N and FU Ori S of 
Am j = 5.6 mag and Am# =4.9 mag that are signif- 
icantly larger than the values estimated by Reipurth & 
Aspin, respectively 4.5 mag and 4.3 mag. 

We thus closely inspe cted the imaging data from 
Reipurth & Aspin (2004). These images were obtained 
using the Infrared Camera for Surveys at the Subaru 
Telescope and are publicly available on the SMOKA 
archival system [^] They consist of a series of images in 
which the brightest star is saturated, and a series of pho- 

10 http://smoka.nao.ac.jp/ 
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tometric calibration frames where the brightest star is at- 
tenuated by a 1/100 neutral density filter. The left panel 
of Fig. [6] illustrates the raw J band IRCS data; the FU 
Ori S component is detected at SNR ~ 3. It is heavily 
embedded in the residual adaptive optics h alo. We can 
reproduce th e delta magnitudes reported in Reipu rth fc| 



Raw PSF with AO halo 



Aspin (2004) by carrying out a gaussian photometric fit 
on these images. However, in an effort to further separate 
the photometric contributions of the AO halo and FU Ori 
S itself, we proceeded to subtract each exposure with a 
centro-symmetric image of itself (right pan el of Fig. [6j 
similar to Fig. 1 in |Reipurth fc Aspin| ( |2004| )) before car- 
rying out the photometric estimation. This process led 
to substantially larger estimates of the magnitude differ- 
ence, Amj = 5.3 mag and Ara# = 5 mag. These values 
are better matched to our PI 640 observations, even con- 
sistent within generous error bars. We then carried out 
the same procedure for the IRCS K' and L' band images 
of FU Ori S. Since the AO halo decreases with wave- 
length, these long wavelength images exhibit SNR > 10 
and thus we found that both methods outlined in Fig. [6] 
yield estimates con s istent with the values published by 
|Reipurth fc Aspm] ([2004) . The results of our photo- 
metric analysis are su mmarized in Table [2] N ote that 
the main conclusion of Reipurth & Aspin (2004), namely 
that FU Ori S is a young star very likely to be associated 
with FU Ori N, relied on red J-H colors. It is thus not 
impacted by our significantly revised photometric esti- 
mates which indicate an even redder object. 



Beck & Aspin ( |2012|) recently reported high spectral 
resolution J H and K band (R ~ 3000) spectra of FU 
Ori S, obtained with the Gemini NIFS integral field spec- 
trograph. They reported a positive SED slope in J and 
a negative slope in H, while our estimates, Fig. [4j yield 
positive slopes in both bands. Although this could be the 
result of extreme near infrared color variability of FU Ori 
S, we have argued against such varaibility for FU Ori N in 
section 2.2, and indeed there may be a technical ex plana- 
tion for the diff erence in spectral slope between the Beck 



|fc As pin (2012) and our observations. While P1640 and 
NIFS are both integral field spectrographs, their design 
tradeoffs are radically different. NIFS prioritizes spec- 
tral resolution (R ~ 3000) over spatial resolution (~ 10 
angular resolution units per spaxel at the shortest wave- 
length). P1640 was designed specifically to explore an 
orthogonal parameter space, with a fine spatial resolu- 
tion in order to mitigate speckle contamination (R ~ 45, 
> 2 spaxels per unit of angular resolution at the shortest 
wavelength). 

In order to investigate the H band slope discrepancy, 
we carried out aperture photometry estimates on our 
P1640 data without any speckle suppression. This led to 
an SED estimate similar to the synthetic s pectrum shown 
Fig. [ 3} which mimics the SED profile in Beck fc Aspin| 
(20121). We thus conclude that the discrepancy in H band 
slope could be due to mid-spatial frequency speckle con- 
tamination in NIFS, that can not be calibrated because 
of the large angular extent of the NIFS spaxels. Alter- 
natively it could be the result of a product near infrared 
variability of FU Ori S and require further monitoring 
of this system. Since the characterization of FU Ori S 



m 



Beck & Aspin (2012) is based on line diagnostics that 



o 




Fig. 6. — In an effort to reconcile the P1640 photomet ry with the 
IRCS data we reprocessed the 2003 epoch published in Reipurth & 
Aspin (2004). Left: raw IRSC PSF. FU Ori S is embedded within 
the adaptive optics halo. Right: Self-calibrated AO PSF, obtained 
by self subtracting each frame with its the centro-symmetric image. 
In the raw image the halo contaminates the photometry of FU 
Ori S and our gaussian fitting analysis yield estimates similar to 
Reipurth Sz Aspin (2004). The same analysis on the Self-calibrated 
AU Jr'S^', where trie halo contamination has been mitigated, yields 
photometric estimates consistent with values obtained on P1640 
data. 

sions. 

4.2. Assembled Spectral Energy Distribution of FU Ori 

S 

We complemented our J and H SED derived from 
P1640 data with published measurements spanning 
0.8 jim to 10 jim. The 0.8 jam point was obtained using 
the delta magnitude of Amp. 8 = 3.96 ± 0.28. r eported 
from speckle imaging by |Kravtsova et al.| (|2008| ) . It was 
placed within the overall SED ol FU Ori S using an ex- 
tr apolation to 0.8 \im of the SE P of FU Ori N published 



Hartmann fc Keny onl (|1996[). We used the K band 



spectrum published by |Keipurth fc Aspin (2004), scaled 
to match the estimated K photometry obtained using our 
analysis of the I RCS data. Comparison of the Reipurth 
& Aspin (2004) and the Beck & Aspin spectral slopes 
through the K band shows that they are similar, if not 
identical, outside of the CO bandhead region. 

The IRCS L band photometry we derived was also 
incorporated in our SED reconstruction. Finally the 
10.2 micron N band point was obtained using a non- 
redundant aperture configuratio n of the Keck segmen t 



Monnier et al. 



;2009) 



take advantage of the high spectral resolution, revising 
their estimated SED slope would not alter their conclu- 



tilting experiment, as reported in 
but with special attention in the data reduction lor our 
purposes to the presence of two point sources. FU OriN 
and FU OriS are separated by 475 ±10 mas and have 
a position angle of 163 ± 2 deg; the N-band flux ratio 
is ~ 10 ± 1.5. Neither component of the binary is itself 
spatially resolved in the 10 fim data. 

To the best of our knowledge there is no other ultravi- 
olet, optical, infrared, or sub- mm data that spatially re- 
solves the pair and therefore could be incorporated into 
our analysis. 

5. ANALYSIS OF THE SPECTRUM AND 
SPECTRAL ENERGY DISTRIBUTION OF FU 
ORIONIS S 

Considering first the spectral line information available 
in our R ~ 45 P1640 data, the J- and H-band spectrum 
lacks any absorption features. This constrains the tem- 
perature to be warmer than ~4000 K since at cooler tem- 
peratures molecular absorption from CO and H2O would 
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T bb = 480 K ; 

A v = 12, 11, h, 9 mag 



; X[/im] 



Fig. 7. — Overall Spectral Energy Distribution of the FU Ono- 
nis system with best fit models. The N-band point is well above 
any conceivable reddened stellar photosphere, and must there- 
fore be due to circumstellar dust around FU Ori S. The colored 
curves are NextGen model atmospheres at the representative ef- 
fective temperature and extinction pairs T e ff = 6500, Ay = 12 
(blue), T eff = 5500, A v = 11 (cyan), T eff = 4500, A v = 9 
(green), T e ff = 3500, Ay = 9 (magenta) all with surface gravity 
log g = 4.5. The red line is a 450 K blackbody with the same 
extinction as the matched photosphere. The sum of the photo- 
spheric and excess contributions are represented by the black lines, 
which are normalized to the data at the P1640 1.7um point. The 
L-band and N-band points can be fit adequately by adding a single 
reddened blackbody of temperature ~ 480 K to reddened photo- 
spheres. 



become apparent even at the low spectral resolution of 
P1640 and in the presence of high extinction. Meanwhile, 
at warmer temperature narrow atomic features may be 
present, but would not be visible at the P1640 resolu- 
tion. The published R ~ 1500 K-band spectrum from 
|Reipurth fc Aspin| ([2004) hints at the presence of Nal at 
2.206 um and per naps Cal at 2.26 /im that the authors 
invoke to advocate for a temperature corresponding to a 
G or K spectral type. This is consistent with their oth- 
erwise featureless K-band spectrum that implies a tem- 
perature warmer than ~4500 K based on the absence of 
CO absorption, and cooler than ~7000 K due to absence 
of BrG absorption (though an actively accreting young 
star could have both CO and BrG absorption filled in 
to the continum level without these lines appearing ex- 
plicitly in emission). Higher signal-to- noise and higher 
sp ectral resolution dat a K-band data recently published 
by |Beck fc Aspjn[ ( [20121 ) more clearly exhibit the photo- 
spheric absorption features in FU Ori S, and suggest a 
spectral type of K5. Notably, these authors find both the 
CO bandhead and BrG in emission, whereas the earlier 



Reipurth & Aspin ( 2004[ ) spectrum showed no emission 
in these regions. It is somewhat unusual for a young 
stellar object to have both clearly present photospheric 
absorption features and strong CO emission. 

The overall SED of FU Ori S from R- to J- to H-band 
is quite red, see Fig. [7) Considering, in addition, the K-, 
L-, and N-bands, the SED is also appreciated as broader 
than a photospheric spectrum. An analysis of the avail- 
able SED was carried out using both blackbody fitting 



and comparison to stellar atmospheres in order to derive 
rough temperature and extinction estimates. We find 
that the N-band point is well above any conceivable red- 
dened stellar photosphere, and must therefore be due to 
circumstellar dust around FU Ori S. The L-band point 
is also likely in excess of the stellar photosphere though 
for very cool, highly reddened atmospheres, an L-band 
excess is not required. Without a more fully populated 
spectral energy distribution the disk can not be charac- 
terized in terms of mass or size, but only in terms of the 
hottest dust closest to the photosphere. The L-band and 
N-band points can be fit adequately by adding to the 
reddened photosphere described below a single reddened 
blackbody of temperature 450 — 480 K . The de-redened 
luminosity of FU Ori N is dominated by its hot inne r 
disk and has been estimated at 226L© ( |Zhu et aL]|2007[ ). 
When combining effective temperature estimates based 
on our reconstructed SED (see below) with a stellar ra- 
dius of ~ li?0 we find that de-reddened luminosity of 
FU Ori N is ~ O.6L . As a consequence, assuming a dis- 
tance of 225 AU between the two components, we find 
that for separations < 10 AU the dusty material around 
FU Ori S is mostly heated by FU Ori S and that further 
away heating from FU Ori N becomes dominant. As the 
temperature from this dust, as per our SED estimate, is 
~ 480K, then it is mostly likely lying within the inner 
AU of the southern component and is most likely heated 
by FU Ori S. 

Moving blueward, the origin of the K-band flux is more 
ambiguous. The shape of both published K-band spectra 
(excluding t he CO bandhead re g ion w hich has different 
slope in the |Reipurth fc Aspin (2004) vs. the [Beck 



Aspin| ( 2012| ) spectra, perhaps indicating the onset of 



CO emission) is consistent with an un-reddened black- 
body of roughly 2200-2500K. The shape of the JH region 
P1640 data, however, requires a much more extincted 
spectrum, for example Av=7 mag for the same 2500 K 
blackbody. This model does not then fit the published 
K-band spectral shape though it does pass through the 
mean K-band flux. Considering alternate blackbodies, 
Av=10 mag is required for a ~5000K temperature black- 
body to fit the JH spectral region, but again this does not 
fit the K-band spectral shape very well, though does still 
match the mean K-band flux. Model atmosphere fitting 
in the JHK wavelength regime introduces sensitivity to 
spectral signatures that would arise in co ol atmospheres. 
For t his purpose we used the NextGen ( |Allard et al. | 
2001[ ) models released in 200^] As seen on Fig.[8j there 
is a general trade-off between extinction and tempera- 
ture through the PI 640 J-band and H-band region, with 
T = 6500,5500,4500, and 3500K models best fit with 
Ay = 12,11,9 and 9 mag, respectively, as illustrated. 
Cooler temperatures exhibit broad molecular absorption 
features throughout the J and H bands, which as noted 
above are not seen. The relatively featureless K-band 
spectrum additionally rules out the 3500K model for 
the same reason that as even cooler models are ruled 
out from the J-and H-band data alone. The ability of a 
model with given temperature and extinction to match 
both the steep positive slope through the J- and H-band 
region and the negative slope through the K-band, de- 
pends on the wavelength at which the model is normal- 

1 1 http: / / phoenix.ens-lyon.fr / Grids/Next Gen / SPECTRA / 



Spectral characterization of FU Orionis S 



9 




log \[fu,m] 




Fig. 8. — Zoom on the region of the SED probed by P1640. Same 
colors as|7] Top: all the models are scaled to match the P1640 spec- 
trum at 1~65 |Um. Bottom: all the models are scaled to match the 
P1640 spectrum at 1.4 /im. Agreement between data and models 
depends on the wavelength at which the models are normalized. 
We find the best-fit for an atmosphere with T = 4000 K and 
Ay = 8, normalized at either the short end of H or the long end 
of J. 

ized to the data. For example, a normalization at the 
far red end of the H-band permits a wider range of tem- 
perature plus extinction models to pass within the er- 
ror bars through the R-band photometry and the J- and 
H-band spectrophotometry, as well as to come close to 
matching the K-band slope, than can be fit to the same 
X 2 -based confidence level when normalizing at shorter 
wavelengths. Shorter wavelength H-band or even J-band 
normalization allows the red end of the H-band as well as 
the K-band to serve as a discriminant between models, 
see Fig.[U Overall, we find the best-fit for an atmosphere 
with T = 4000K and Ay = 8, normalized at either the 
short end of H or the lon g end of J. This effecti v e temper- 
ature is consist ent with |Beck fc Aspin (2012); Reipurth 
fc Aspin (2004). our PI 640 observations provide tighter 
constraints on this quantity along with an estimate of 
the extinction in the line of sight of FU Ori S. 

6. NEW RESULTS ON FU ORI S IN THE 
CONTEXT OF FU ORI N 



In this paper we have established that FU Ori S is 
physically related to FU Ori N and the pair forms a true 
binary system. The P1640 astrometry over 2 epochs 
combined with literature information establishes com- 
mon proper motion of the pair. Our spatially resolved 
spectrophotometry in JH in combination with the lim- 
ited amount of data at other wavelengths that spatially 
resolves the binary enables for the first time an assess- 
ment of the N/S pair of young stars in the FU Orionis 
system. The northern component is a long recognized FU 
Ori star that is visible at optical wavelengths and well- 
studied in the context of the accretion disk-dominated 
paradigm for FU Ori objects, the class for which it is 
the prototype. Previous modelling work has estimated 
that the underlying star is a 0.3 to 0.5 M star behind 
Ay = 1.5 — 2.4 mag. The southern component, which we 
have characterized here, appears to be the more heavily 
reddened (Ay = 8 — 12) component of the system, and 
perhaps even the slightly hotter (~ 4000-6500 K) and 
therefore more massive (> O.5M ) component, depend- 
ing on the adopted temperature. Both components have 
infrared excess though there are no spatially resolved flux 
measurements beyond those at 10.7 jam that are newly 
reported here. 

Notable is the much higher extinction required in or- 
der to explain the optical to near-infrared spectral en- 
ergy distribution of FU Ori S, compared to the more 
moderate reddening inferred from other studies towards 
FU Ori N. This difference should be put in perspective 
with extinction estimates obtained using X-ray observa- 
tions of this system. Two separate line of sight absorp- 
tions are needed to fit the FU Ori X ray spectrum to 
thermal emission m odels. Building upon earlier work 
with XMM-Newton dSkinner et 



Idmg upo n earlier work 
al. |2006D , |Skinner et ah] 

ita to show that the FU 



d2010| ) used Chandra AC1S-S data 
Ori spectrum is fit well with (1) a highly absorbed hot 
plasma component (kT ~ 4.5 keV; N H — 10 23 cm -2 ), 
that is spatially coincident with the primary, and (2) a 
cooler, less absorbed thermal plasma (kT ~ 0.5 keV; 
NH ~ 10 22 cm -2 ) which shows a small centroid offset 
towards FU Ori S. They suggest that the latter soft com- 
ponent is most likely due to FU Ori S, which contributes 
at least to 50% of the soft X-ray counts. Table [3] re- 
ports the values of th e estimated optica l extinction for 
the FU Ori N-S pair, [Zhu et al.| fl2007| ) and this work 
respectively, converts them to hydrogen column density, 
and compares these indirec t Nh estimates with the di- 
rect X-ray measurement in | Skinner et al. (2010). Esti- 
mated hydrogen column densities in the direction of FU 
Ori S obtained using the two methods are of the same 
order of magnitude and consistent within generous error 
bars. The independent estimate derived from our obser- 
vations thus reinforces the notion of the presence of a 
cool plasma embedding FU Ori S. Moreover it could po- 
tentially be useful to externally constrain the soft X-ray 
compo nent in the global fit presented in | Skinner et al.| 
(2010). Such an exercise would yield to fitting the Chan- 
dra counts with a large hydrogen column density for the 
h ard X-ray counts, eve n larger than the values presented 



Skinner et al. (2010). This would further reveal a dis- 



crepancy between optical extinction measurements and 
X-Ray absorption the direction of FU Ori N even larger 
than the factor > 10 exhibited in Table (3) This excess 
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X-ray absorption indicates a large gas-to-dust ratio that 
could be explained by non-uniformities in the geometry 
of FU Ori N' s disk flKravtsova et al.||2007||Kley and Lin 
|1996| > or the presence ol winds. A full understanding ol 



this paradigm will require further researc h, and we re- 
fer the reader to the detailed discussion in Ski nner et al.l 
(2010) for a thorough discussion of potential avenues to 
identify such phenomena. 

The infra red excess, established he re and implied by 
the data in |Reipurth fc Asprn] (|2004|) as well as the CO 
plus BrG emission reported by [beck & Aspin| ([2012]) 
(though m ysteriously not present in the |Keipurth fc As- 
pin (2004) spectrum), strongly suggests that the FU Ori 
S is also a young star/disk system. However, the lack 
of cool molecular features and the implied temperature 
of at least 4000 and perhaps more than 5000 K actu- 
ally render FU Ori S the more massive object, in other 
words, the primary of the system. The observed opti- 
cal and infra-red brightness ratios can then be explained 
by both the excess luminosity of FU Ori N due to its 
outburst state since the 1936 eruption, and the high ex- 
tinction in the line of sight towards FU Ori S. Dust in 
the hot inner disk around FU Ori N cannot account for 
this extinction since the pair separation of ~ 225 AU is 
large r than published estimates of the outer radius of this 
disk (Zh u et al.|2007 ). Perhaps, however, the pair is em- 
bedded in the same circumbinary envelope with higher 
extinction in either a clumpy envelope or the foreground 
cloud present along the line of sight to the S component 
than towards the N component. Indeed, the pair could 
be a binary in the process of fragmentation. The source 
may thus be analogous to the Z CMa binary system, 
where the optically visible source is also a candidate FU 
Ori star and the embedded companion is the more mas- 
sive object, called a Herbig Ae/Be star in the literature. 
The Z CMa system is disc ussed in a separate pa per from 

t he P1640 collaboration (([Hinkley et al. |2Q12|» . 

Reipurth fc Aspin|(|2^ fc Syer ( HPS) ; |Clarke 



et al. ( p005[ ) and |Bonnell fc Ba stien ( 1902]) have proposed 



binary interactions as a suitable trigger for the FU Ori 
outbursts, causing gravitational instabilities that raise 
accretion rates by 3-4 orders of magnitude over normal 
T Tauri quiescent accretion levels. Problem s with this 
origin of th e outbursts (see |Pfalzner (2008); [Forgan fc 



Rice (2010)) include the short predicted rise times, the 



higher order multiplicity suggested as necessary for aperi- 
odic instabilities than the mere binarity that is observed, 
and finally that the gravitational instability would affect 
the entire disk, not just in the inner fraction of an AU 
that is thought to be participating in the outburst. Note 
that in our case the separation between FU Ori N-S is 
too wide for FU Ori S to be responsible for the 1935 
outburst. 

7. CONCLUSION 

We have presented near-infrared Integral Field Spec- 
trograph observations of FU Orionis: our data spatially 
resolves the FU Ori N and FU Ori S binary components 
throughout the J- and H-bands, at R~45. Our astromet- 
ric analysis unambiguously establishes common proper 
motion FU Ori N and S form a binary pair. Our Spectral 
Energy Distribution analysis suggests that the southern 
faintest source in this system might actually be the more 
massive component. 



Our observations allowed us to retrieve the SED of FU 
Ori S, which was thus far poorly constrained. In order 
to carry out unbiased spectro-photometric estimates in 
the pres ence of speckles we applied the damped LOCI al- 
gorithm |Pueyo~eFffl ([2011]), a reduction method specifi- 
cally designed lor high contrast science with IFS. This is 
the first communication reporting the high accuracy, on 
a faint astronomical source, of this technique pioneered 
by the Project 1640 team. 

We combined the PI 640 data with all other spatially 
resolved information (R-band, K-band spectroscopy, L- 
band, and N-band) available from the literature or from 
re-analysis of previously existing data. L and N-band 
point are well above any conceivable reddened stellar 
photosphere, and must therefore be due to circumstellar 
dust around FU Ori S, whose average temperature we es- 
timated at ~480 K. The spectral energy distribution of 
FU Ori S is very red from optical wavelengths through 
the H-band with a turnover thereafter such that the Ri- 
band spectral slope is blue. This combined data set is 
best explained by an underlying star in the 4000 — 5000 
K temperature range behind 8 — 12 mag of visual extinc- 
tion. 

This value of the extinction is in good agreement 
with estimates obtained using previously published X-ray 
measurements. In particular our estimate is consistent 
with the hydrogen column density derived using Chan- 
dra soft X-r ay counts, which have been mostly attributed 
to FU Ori S |Skinner et al.| ( |2010| ). Our independent con- 
straint on tne extinction in the direction of FU Ori S 
helps in turn to further confirm the anomalously large 
excess hard X-ray absorption towards FU Ori N. 

The high visual extinction in the line of sight of FU 
Ori S, most likely due to the geometry of a potential cir- 
cumbinary dust shell, and is moreover responsible for the 
relatively high contrast ratio with respect to FU Ori N. 
Our SED analysis allows us to combine estimates and of 
visual extinction and effective temperature, and yields 
values of the latter twice larger than previously pub- 
lished. FU Ori S is thus much more massive than pre- 
viously believed and, with a lower bound for the mass 
estimate of > 0.5Mq, it is probably the most massive 
component of this system. This source may thus be anal- 
ogous to the Z CMa binary system. A more detailed 
comparative analysis of these two objects, putting in per- 
spective the quiescent and outburst near infrared SEDs 
of each components in both FU Orionis and Z CMa, is 
a promising avenue to better unravel role of binarity in 
eruptive protostars. 
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TABLE 1 

RELATIVE ASTROMETRY OF THE FU ORIONIS SYSTEM 



Epoch ARA ADEC Referen ce 



Oct 27 th 2002 


0.16" 


±0.03" 


-0.47" 


± 0.03" 


Wang et al. (2004) 


Dec 15 th 2003 


0.147" 


± 0.006" 


-0.470" 


± 0.006" 


Reipurth & Aspin (2004) 


Feb 20 th 2005 


0.139" 


±0.017" 


-0.454" 


± 0.017" 


Monnier et al. ( 2009 ) 


Jan 23 rd 2008 


0.137" 


± 0.008" 


-0.474" 


± 0.008" 


Kravtsova et al. ( 2008 ) 


Mar 17 th 2009 


0.158" 


± 0.009" 


-0.465" 


± 0.009" 


This work 


Dec 12 th 2011 


0.15" 


± 0.02" 


-0.45" 


± 0.02" 


This work 



TABLE 2 

PHOTOMETRY OF THE FU ORIONIS SYSTEM 



Filter 


A c 


FU Ori N 


FU Ori S 


Instrument 


J 


1.250 


6.519 ±0.015 


11.98 ±0.19 


P1640^ 


J 


1.250 


6.519 ±0.015 


11.69 ±0.22 


IRCS W 


H 


1.635 


5.699 ±0.029 


10.36 ±0.17 


P1640 ^ 


H 


1.635 


5.699 ±0.029 


10.14 ±0.19 


IRCS W 


K' 


2.121 


5.259 ±0.023 


9.35 ±0.15 


IRCS W 


L' 


3.770 


4.180 ±0.039 


8.09 ±0.16 


IRCS W 


N 


10.7 


2.2 ±0.02 


4.7 ±0.13 


Keck Segment masking ( c ) 



Note: (a) Band averaged phot ometry based on the P1 640 spectra. 

(b) IRCS program in |Reipurth fe Aspin| j2004| reprocessed 
in this work. 

(c) Data from Monn ier et al.| ( |2009| ) reprocessed in this 
work. 



TABLE 3 

EXTINCTION IN THE FU ORIONIS BINARY PAIR 



Line of sight 


Ay in mag ^ 


N H = f(A v ) in 10 22 cm" 2 W 


N H in 10 22 cm ~ 2 ^ 


FU Ori S 


10 ±2 


2 ±0.4 


1.1 ±0.5 


FU Ori N 


2 ±0.5 


0.39 ±0.11 


10 ±5 



Note: (a) |Zhu et al] ( [2010| for FU Ori N, this work for FU Ori S. 

(b) Uon version carried out using \fty :| 

//_ ftp.astro.princeton.edu / draine / dust / mix I kextglbedow D m W3 .1q0 

(c) Be st ht tor the two components model in | Skinner et al.] 
|2010| ). 
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